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I .  INTRODUCTION 


The  question  of  phase  retrieval,  that  is,  of  the  determination  of 
the  phase  of  an  incoming  optical  signal,  is  an  extremely  important 
problem,  one  which  can  be  approached  in  several  different  ways. 

One  of  the  most  practical  methods  of  phase  retrieval  is  to 
attempt  to  deduce  the  phase  of  the  wavefield  from  a  knowledge  of  the 
intensity  of  the  field  in  two  different  planes:  the  aperture  plane 
and  focal  plane  of  a  thin  lens. 

In  this  report  we  will  investigate  several  aspects  of  the  unique¬ 
ness  of  the  phase  retrieval  performed  by  this  method,  that  is,  whether 
or  not  the  method  allows  one  to  uniquely  determine  the  phase  of  the 
incoming  wavefield. 

The  method  of  phase  retrieval  under  consideration  can  be  described 
more  precisely  in  the  following  way.  Let  the  field  in  the  aperture 
plane  of  the  lens  (which  we  take  to  be  the  plane  z  =  0)  be  written  as 

V(P0.t)  -  Re  [Ba(J0>  e  lut]  , 

-Re  [|0,<?0)|  ,  (1.1) 

where  =  (xg»yg)  is  the  two-dimensional  position  vector  in  the 
aperture  plane,  U  (pn)  is  the  complex  amplitude  of  the  field  in  that 
plane,  w  is  the  temporal  frequency  of  the  field  and  Re  denotes  that 
the  real  part  is  to  be  taken.  The  method  in  question  strives  to  de¬ 
termine  the  phase  kW(pg)  from  the  knowledge  of  the  aperture  plane 
intensity,  |Ua(Pg)|  and  the  focal  plane  intensity,  |Uf(p)|^.  p  =  (x,y 
is  the  position  vector  in  the  focal  plane.  From  here  on  ju|  will  be 


referred  to  as  the  amplitude  of  the  field. 

Numerical  experiments  which  successfully  retrieve  the  phase  from 
aperture  plane  and  focal  plane  intensities  have  been  performed.  Gon¬ 
salves^  treated  the  one-dimensional  focusing  problem.  He  found, 
for  several  test  cases,  that  his  algorithm  always  converged  to  either 

the  correct  phase  or  its  "twin"  (the  twin  solution  will  be  discussed 

2 

in  Section  V) .  Southwell  treated  the  two-dimensional  problem  for  the 
case  of  uniform  aperture  plane  intensity,  |u  ( p.)|  ^  =  1,  and  a  phase 
described  by  means  of  the  nine  Zernike  polynomials  corresponding  to 
tilt,  defocus,  astigmatism  and  coma.  He  found,  for  several  test  cases, 
that  his  algorithm  produced  the  correct  phase. 

Neither  of  these  papers  made  a  definitive  statement  about  the 
general  problem  of  uniqueness,  that  is,  about  whether  or  not  there 
is  only  one  possible  phase  function  consistent  with  the  two  known  in¬ 
tensities.  Until  this  point  is  clarified,  one  cannot  be  confident  that 
the  phase  one  calculates  is  the  one  occuring  physically.  In  the  one¬ 
dimensional  problem,  some  interesting  results  on  the  question  of 

3  A  5  6 

uniqueness  have  been  obtained.  ’  ’  ’  However,  the  two-dimensional 
problem  is  not  a  straight-forward  extension,  because  physical  wave- 
functions  cannot,  in  general,  be  subjected  to  the  method  of  separation 
of  variables  and  furthermore  lenses  have  circular  apertures,  not 
square  ones.  Indeed,  recently  the  question  of  uniqueness  for  the 

7  8 

problem  we  have  described  above  has  been  debated  in  the  literature.  ’ 

More  recently  two  proofs  of  uniqueness  have  been  obtained.  Foley 
9 

and  Butts  have  shown  that  if  the  wavefield  intensity  across  the  lens 
is  uniform  and  the  phase  is  represented  by  means  of  the  nine  lowest 
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order  Zernike  polynomials,  that  the  phase  can  be  determined  uniquely 
(except  for  the  well-known  twin  solution  ambiguity)  from  its  aperture 
plane  and  focal  plane  intensities.  Foley ^  has  extended  this  proof  to 
include  aperture  plane  intensities  which  are  Gaussian. 

The  purpose  of  this  report  is  to  investigate  two  aspects  of  the 
uniqueness  question.  First,  using  a  specific  wavefield  (an  astigmatic 
Gaussian  beam),  we  will  investigate  the  role  that  diffraction  plays  in 
the  question  of  uniqueness.  (By  "diffraction"  we  mean  that  due  to  the 
finite  size  of  the  lens,  not  diffraction  due  to  unobstructed  propagation 
in  free  space).  It  will  be  shown  that  this  diffraction  is  responsible 
for  the  uniqueness  which  arises.  Secondly,  we  will  discuss  the  twin 
solution  ambiguity  and  propose  a  method  for  eliminating  it. 


II.  PHASE  RETRIEVAL  FROM  INTENSITY  MEASUREMENTS 


In  this  section  the  particular  method  of  phase  retrieval  which  we 
are  considering,  namely  the  determination  of  the  phase  of  an  incoming 
complex  amplitude  wavefield  from  the  knowledge  of  the  intensities  in 
the  aperture  plane  and  focal  plane  of  a  thin  lens,  will  be  described 
in  a  more  precise  mathematical  fashion. 

2 . 1  Precise  Mathematical  Statement  of  the  Uniqueness  Problem 
Let  the  wavefield 

V±(r , t)  «  Re  {Ui(?)e'iu)t}  ,  (2.1) 

where  r  =  (x,y,z),  be  traveling  towards  positive  values  of  z.  (From 
here  on  the  e  il0t:  will  be  omitted,  since  the  field  at  any  position  in 
space  will  have  this  form  of  time  dependence.  We  will  now  deal  with 
complex  amplitudes  only  and  refer  to  them  as  the  fields.)  Furthermore 
let  this  field  be  incident  upon  a  thin  lens  of  focal  length  f  and 
radius  L,  which  is  located  at  z  =  0  and  surrounded  by  a  stop.  In  this 
plane  the  incident  field  can  be  written  in  terms  of  its  modulus  and 
phase  as 

U^Pq.O)  =  |Ui(?0,0))eikW(p0),  (2.2) 

where  =  two  dimensional  position  vector  in  the  plane  z  =  0,  k  =  w/c 
and  W(p  )  is  a  real  function  of  p^.  See  Figure  2.1.  The  aperture  plane 
field  (which  includes  the  effect  of  the  stop)  is 

Ua(P0)  =  ^(Pq.O)  circ  (p0/L)  ,  (2.3) 


where  pQ  =  |pQ|  and 


1,  P0  1  L 

circ  (Pq/L)  =  .  (2.4) 

0,  P0  >  L 

The  field,  U^(p) ,  in  the  focal  plane  can  be  found  by  using 
diffraction  theory.  It  turns  out  to  be  proportional  to  the  Fourier 
transform  of  the  aperture  plane  field^, 

ikf  ikp2/2f 

Uf(p)  =  - - -  l?a(p/Af)  ,  (2.5) 


where  A  is  the  wavelength  of  the  light  and 


-v  ->  *  -2iti(p/Af)*p  , 

lUp/Af)  =  _//U  (pn)  e  °dP„. 


(2.6) 


It  follows  from  Eq.  (2.5)  that  the  intensity  distribution  in  the  focal 
plane  is 


If(p)  =  tuf (P) 1 2  . 

4?  }U  (p/Af)  | 2  .  (2.7) 

aV  a 

The  phase  retrieval  problem  can  be  stated  as:  given  the  intensity 
distribution  in  the  focal  plane  and  aperture  plane  of  a  lens,  can  we 
then  uniquely  determine  the  phase  of  the  aperture  plane  field?  Or  in 
other  words,  is  it  possible  for  two  (or  more)  aperture  plane  fields  with 
equal  aperture  plane  intensities  but  different  phases,  e.g., 
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ikW  (p  ) 

UaB<pO>  =  *  ’ 

to  produce  equal  focal  plane  intensities?  The  equality  of  the 
aperture  plane  intensities  demands  that 


^aA^i2  =  ^aB^O*1  2  '  f°r  a11  P0  * 


(2.8) 


The  equality  of  the  focal  plane  intensities  demands  that 


[UaA(pAf)!2  -  |Ua*(p/Af)|  2  ,  for  all  p, 


(2.9) 


2.2  The  Twin  Solution 

13  12 

It  is  well-known  *  ’  that  if  the  aperture  plane  amplitude  is 
symmetric  with  respect  to  inversion,  i.e.. 


IV-?o> 


lua<p0>l  »  for  a11  po’ 


(2.10) 


then  there  are  at  least  two  fields  which  will  have  the  same  intensities 
in  both  the  aperture  plane  and  focal  plane.  They  are  the  field  itself, 

ikWA(Po) 


UaA<P0)  =  KM  6 


(2.11) 


and  its  twin  , 


\*$0>  •  l“a<»0) 


-ikWA(-Po) 


(2.12) 


These  two  fields  have  the  same  aperture  plane  intensities  by  defi¬ 
nition.  Furthermore,  as  we  will  shown  in  Section  V, 


Vp/Af)  =  [^(p/Af)]*  , 


(2.13) 


L. 


i 


8 


hence,  from  Eq.  (2.7),  we  see  that  they  have  the  same  focal  plane  in¬ 
tensities  as  well.  In  Section  V  we  will  discuss  the  twin  solution 
further.  In  Sections  III  and  IV  we  must  keep  in  mind  that,  since  we 
will  be  dealing  with  fields  which  obey  Eq.  (2.10),  the  twin  solution 
is  always  a  possibility. 
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III.  THE  ROLE  OF  DIFFRACTION  BY  LENS  APERTURE  IN  THE  QUESTION  OF 

UNIQUENESS.  PART  I:  RADIUS  OF  BEAM  RADIUS  OF  LENS 
8 

Devaney  in  his  paper  On  the  uniqueness  question  in  the  problem 

of  phase  retrieval  from  intensity  measurements"  suggested  that  although 

the  phase  retrieval  problem  that  we  are  discussing  was  shown  by  Robinson^ 

to  be  nonunique  in  the  limit  of  geometric  optics,  it  does  however  appear 

to  possess  a  unique  solution  (except  for  the  twin  solution)  within  the 

framework  of  diffraction  theory.  Indeed,  uniqueness  (except  for  the 

twin  solution)  for  the  problem  we  are  considering  has  been  proven  by 
9 

Foley  and  Butts  . 

In  this  and  the  next  section  we  will  investigate  the  role  of  the 
diffraction  caused  by  the  lens  aperture  in  the  question  of  uniqueness. 

We  will  study  a  particular  type  of  incident  field,  namely  an  astigmatic 
Gaussian  beam.  In  this  section  we  will  assume  that  L,  the  radius  of 
the  lens,  is  much  larger  than  the  radius  of  the  incident  Gaussian  beam, 
so  that  there  is  no  diffraction  due  to  the  lens  aperture.  We  will  find 
that  this  leads  to  a  great  deal  of  nonuniqueness,  i.e.,  several 
beams  with  different  phases  can  have  the  same  intensity  in  hoth  the 
aperture  and  focal  planes.  In  the  next  section  we  will  treat  the  case 
where  the  beam  radius  is  not  negligible  compared  to  L  and  we  will  see 
that,  except  for  the  twin  solution,  all  the  nonuniqueness  is  removed. 


natic  Gaussian  Beams 


Let  us  now  investigate  the  case  where  the  incident  fields  is  an 


astigmatic  Gaussian  beam.  The  incident  field  can  be  written  as 


-  .  ikw(p0> 

Vv°)  *  e 

(3.1) 

where 

.  .+  ,  [2  1  -p0  ^0 

A(poJ  JiT  wQ  e 

(3.2) 

is  a  normalized  Gaussian  amplitude,  and  the  phase  function  W(Pq)  can 

be  represented  by  the  expansion, 

5 

W(pn)  -  E  a  Z  (p  /L)  . 

j*=l  J  J 

(3.3) 

The  aj's  are  expansion  coefficients  and  the  Z^'s  are  a 
13 

Zernike  polynomials 

set  of  modified 

z1(pq/l)  »  xq/L  , 

(3.4a) 

Z2(pQ/L)  =  yQ/L  , 

(3.4b) 

Z3(pq/L)  =  2(p 0/L)2  -  1  , 

(3.4c) 

z4(p0/l)  -  (x0/l)2  -  (y0/D2  , 

(3.4d) 

Z5(pq/L)  =  2(x0/L)(y0/L)  . 

(3.4e) 

Each  of  these  terms  has  a  meaningful  physical  interpretation.  The 
first  two  terms  represent,  respectively,  the  tilt  in  the  x  and  y  directions. 
The  third  term  corresponds  to  defocus  and  the  last  two  terms  represent, 
respectively,  0°  and  45°-astigmatism. 

Substituting  Eq.  (3.2),  (3.3),  and  (3.4)  into  Eq.  (3.1),  we  then 


re 


11 


-v  [2  1  -P0  /»o  (  IVo 

ui(p0'0)  e  expjik[— + 


Vo  ,  V2P0  -  » 


a4UQ2-  yQ2)  +  ^5*">' 


*1!  ■ 


(3.5) 


Furthermore  for  the  sake  of  mathematical  convenience  we  can  drop  the 
linear  terms  in  Eq.  (3.5),  since  their  only  effect  would  be  to  shift 
the  focal  plane  field  distribution  by  x  =  -a^f/L  and  y  =  -a^f/L. 

Eq.  (3.5)  can  then  be  written  as 


Vv0) 


-P-L 

J’“o 


exp{ik[a3(pQ2  -  \  L2) 


,  2  2 
a4(x0  ‘  y0 


+  2acx»y_  ]}  ,  (3.6) 


where 


2a3/L2  , 

(3.7a) 

a^/L2  , 

(3.7b) 

a5/L2  . 

(3.7c) 

Throughout  the  rest  of  this  report  our  incident  field  will  be  of 
the  form  of  (3.6).  Note  that  measures  the  defocus,  measures  the 

0°  astigmatism  and  measures  the  45°  astigmatism.  The  tilt  terms 
(a^  and  a^)  are  needed  only  if  the  focal  plane  intensity  pattern  is 
not  centered  about  x  =  y  =  0. 

3.2  Propagation  of  an  Astigmatic  Gaussian  Beam  Through  a  Thin  Lens 
(Infinite  Aperture) 

3.2.1  Astigmatic  Gaussian  Beam  with  ■  0 
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1/R  =  1/R-  -  1/f 

yo  yo 


(3.11b) 


In  Eq.  (3.10)  we  have  used  the  fact  that  ^(Pq)  =  ^or  t*le 

case  of  an  Infinite  aperture.  If  we  now  use  the  complex  radius  of 
curvature  notation^,  Eq.  (3.18)  can  be  rewritten  as 


where 


U(p0,0+)  - 


(3.12) 


1/q  =  1/R+  +  iA/irW  2  , 

x0  X0  0 


(3.13a) 


1/q  =  1/R+  +  iX/nU  2  . 

yo  yo 


(3.13b) 


3.2.2  Beam  propagation  (z  >  0) 

We  are  now  ready  to  calculate  the  field  in  any  plane  z  =  constant 
>0.  We  will  use  the  Fresnel  diffraction  integral, 


ikz  _  „  ik(p-pQ)2/2z  2 


U(p,z)  =jzr-%  u(V0+)e 


d  P0’ 


ikz  ik  2z 

e  e _ 

ilz 


+  ikpQ2/2z  -ikp*pQ/z  2 
f/U(p0,0  )  e  0  e  0  d  pQ.  (3.1- 

It  follows  in  a  straightforward  manner, after  substituting  (3.12) 
into  (3.14),  that  the  resulting  field  is 

■1.1, 2/u  2t,\  iki2",~  -  2 2 

u(p,z)  =J|- 


4(z)  -X2/ W  2(z)  ikx2/2R  (z)  -y2/W  2(z)  iky2/2R  (z) 

T7T) - V~d - (3’15) 

x  y 


where 


ux(2>  ■  "o  /^TT^T- , 


(3.16a) 


VZ)  =  W0  J(1  +  z/Ry  >  +  <z/zR> 

(1  +  z/R+  )2  +  (z/z_)2 

X  K 

R  (Z) - - - - - — 

(1/R  )(1  +  z/R+  )  +  (z/zp 

x0  Xq  k 


(3.16 


(3.17 


(1  +  z/R+  )2  +  (z/z  )2 
y0  K 

R  (z) - - - - - 5“ 

y  (1/R*  )(1  +  z/r: )  +  (z/zf) 


(3.17 


zR  =  irWQ  /X  , 


(3.18 


ij)(z)  *  k(z  -Va  L  )  (z)  (z) 

j  x  y 


(3.19 


*x(z) 


-1  XhV0 


(1/ R,  +  1/z) 
0 


(3.20 


i(/y(z> 


-1  A/*W0 


(1/R  +  1/z) 


(3.20 


The  intensity  is  therefore 


I(p,z)  -  U  (p  » z)  U(p,z)  , 


-2x2/V  2<*>  -2y2/W  2(z) 


“x(‘> 


W  (z) 

y 


(3.21 


From  Eq.  (3,21),  we  see  that  for  any  particular  z,  the  isophotes  (con 
tours  of  equal  intensities)  are  specified  by 


x2/W  2(z)  +  y2/W  2(z)  =  constant. 
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r 


which  are  ellipses  which  can  degenerate  Into  circles  If  W  (z)  =  W  (z) . 

x  y 

Circular  isophotes  (If  they  occur)  are  observed  at  only  one  particular 
value  of  z  >  0.  This  circular  image  is  called  the  circle  of  least 
confusion. 


3.2.3  Propagation  of  a  General  Astigmatic  Gaussian  Beam 

Let  us  now  investigate  the  case  where  is  not  generally  zero. 

18 

(We  will  be  following  formalism  similar  to  that  developed  by  Cook.  ) 
The  field  on  the  left  side  of  the  lens  can  then  be  written  as 


ba-°> -Jh 


2  1  W  lkH<»0> 

-  77-  e  e 


(3.22) 


where 


W(P0)  "  a3(P02"  **  L)2+  o4<x02  -  yQ2>  +  2a5x0y0  *  (3*23) 


As  in  Sec.  3.2.1  we  can  now  write  the  field  on  the  right  hand 
side  of  the  lens  as  2 

-  +  +  -ikP0  /2f 

U(P0.°  )  «  Ua(pQ)  e 


(2  1  ->0IV0  -lk<>02/2f  lkW<?0> 

“Jx  e  6  6 


(3.24) 


Using  the  Fresnel  diffraction  integral,  the  field  for  z  >  0  is 


U(p,z) 


ikz  lkp  /2z 


e  e 
iXz 


jo  +  lkP0  /2z  T-ik  1 

J/U(pQ,0  )  e  exp[— (xxQ  :  yy0)J  dxQdy0 


(3.25) 
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Substitution  of  Eq.  (3.24)  into  (3.25)  yields 

v  eikz  ikp2/22  u  -P02/W02  “ikP02/2f  lkP02/2z  lkW(iV 

U(p,z)  = - jjj - 

xexp  (xxQ  4-  yy0)]  dxQdy0  .  (3.26) 

2ika5X0y0 

The  term  e  makes  this  integral  nonseparable  and  different 

from  the  integral  we  encountered  earlier.  Thus,  we  have  to  go  to  a 
rotated  coordinate  system  in  which  this  term  vanishes  in  order  to  make 
the  integral  separable. 

Let  us  now  rotate  the  aperture  plane  coordinates  and  the  observation 
plane  coordinates  counterclockwise  by  an  angle  0  about  the  z-axis.  In 
order  to  do  this  we  now  define  a  new  aperture  plane  coordinates,  (x^.y^) , 
and  a  new  observation  planes  coordinates,  (x,y) ,  according  to  the  trans¬ 
formation  equations, 

xo  =  xocose  +  Vine 

yQ  «  -xQsin0  +  yQcos0  ,  (3.27) 

and 

x  =  xcos0  +  ysinQ  , 

y  =  -xsin8  +  ycos6  .  (3.28) 

The  old  coordinates  are  then  given  in  terms  of  the  new  coordinates 
according  to 

Xq  =  xQcos0  -  y0sine 

yQ  -  x0sin6  +  y0cos8  ,  (3.29) 
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(3.35a) 
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Uq-  =1/4  +  iX/ttW  2  , 

yo  yo  u 

(3.35b) 

1/R-+  =  2(a-  +  a.)  -  1/f  , 

xo  J  H 

(3.36a) 

1/R_+  =  2(a3  -  a4)  -  1/f  . 

(3.36b) 

Eq.  (3.32)  has 

the  same  form  as  Eq.  (3.14)  with  Eq.  (3.12) 

substituted 

into  it.  It  then  follows  that 

U(x,y,  z)  =  If 

l?(z)  -x2/W_2( z)  ikx2/2R_(z)  -y2/W_2(z)  iky2/2R_(z) 

eexe  xeye  y 

Jv?  iv*7 

(3.37) 

where 

VZ) 

■"o  J(1  +  +  ■ 

(3.38a) 

W_(z) 

•”o  ]<1+  z/E5o)2+  U/Zr>2  ’ 

(3.38b) 

V2> 

(1  +  z/R_V  +  (z/z  )2 
xo  R 

(1  +  z/R_+) /R_+  +  (z/z?) 

xo  xo  R 

(3.39a) 

Bj(z) 

(1  +  z/R_V  +  (z/zp)2 
yo  R 

(1  +  z/R-+)/R-+  +  (z/z2) 

yo  yo  K 

(3.39b) 

ii(z) 

2 

=  k(z  ~  4a3L  )  “  *5^  ~  , 

(3.40) 

and 
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i/>_(z)  =  tan 


t|»—  (z)  =  tan 


(3.41a) 


(3.41b) 


The  intensity  in  the  rotated  frame  is  then  given  by. 


I(x,y,z) 


-2x2/wJ(z)  -2y2/W_2/(z) 

2  e  x  e  y 

v  w^(z)  Wy(z) 


(3.42) 


It  follows  then  that  for  any  particular  z,  the  isophotes  are  generally 
ellipses  (Figure  3.1)  in  a  rotated  plane.  It  also  follows  that  circular 
isophotes  (if  they  occur)  are  observed  at  only  one  particular  value  of 
z  >  0. 


3.3  Nonuniqueness  in  the  Case  of  an  Infinite  Aperture 

The  uniqueness  question  asks:  can  two  (or  more)  astigmatic  beams 
with  equal  aperture  plane  intensities  but  different  phases,  e.g., 


„  <*> . ix ± rp» /u° 

aAv*V  Ji  W. 


-pn2/Wn2  ika3PQ2  ikaA(x02  -  y  2)  ika  2x  y 
e  e  e 


2,..  2 


2  2. 


„  ,  .  p  I  "P0  /W0  lke3P0  ik64(xO  “  y0  }  ik652Vo 

UaB(p0  Jx  ^  6 


generate  identical  focal  plane  intensities?  It  follows  from  Eq.  (3.42) 
that  the  two  fields  will  have  the  same  focal  plane  intensities  if  their 
x  and  y  beam  widths  are  equal  in  the  focal  plane, 
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W_A(f)  =  W— B(f)  ,  (3. A3) 

W_A(f)  =  wJ*(f)  ,  (3 .44) 

and  the  rotation  angle  of  the  focal  plane  isophotes  are  equal, 

6  A  =  6  B  .  (3.45) 


By  using  Eq.  (3.38),  (3.36)  and  (3.34),  Eq.  (3.43)  -  (3.45)  can  be 
written  as 


(3.46) 

(3.47) 


(3.48) 


Tables,  1,  2,  and  3  list  all  the  possible  solutions  to  this  set 
of  equations.  Table  3.1  gives  all  the  possibilities  for  the  case  of 
=  0,  Table  3.2  gives  all  the  possibilities  for  the  case  of  =  0 
and  Table  3.3  gives  all  the  possibilities  for  the  case  where  both 
and  are  not  equal  to  zero. 

It  is  apparent  from  Table  3.1  that  in  this  case  there  are  an 
infinite  number  of  fields  which  will  duplicate  the  intensity  created 


by  the  field  /  0,  =  0,  since  we  are  free  to  choose  6^  to 

2  2 

be  anything,  as  long  as  =  +  la ^  -  8^  .  The  same  is  true  for  some 


of  the  solutions  in  Table  3.2.  This  occurence  of  an  infinite  number 


of  solutions  is  caused  by  the  fact  that  all  the  fields  created  by  the 
solutions  in  Table  3.1  and  Table  3.2  produce  circular  isophotes  (see 
Appendix  A)  in  the  focal  plane.  Hence,  since  the  aperture  plane  field 
is  also  circularly  symmetric  any  rotation  of  this  incident  field  (which 
corresponds  to  a  different  8^  and  therefore  a  different  8^)  will  leave 
the  focal  plane  field  unchanged  as  long  as  8^  and  8^  are  related  in 
the  manner  stated  above.  Thus,  in  practice,  if  you  have  an  astigmatic 
beam  incident  upon  a  thin  lens  where  L>>W^,  and  you  see  circular  iso¬ 
photes  in  the  aperture  plane  and  focal  plane  there  is  no  hope  of 
determining  the  phase  of  the  incident  field. 

Table  3.3  corresponds  to  fields  which  create  elliptical  isophotes 
in  the  focal  plane.  Here  there  are  only  four  possible  solutions:  two 
distinct  solutions,  each  with  its  own  "twin". 

In  order  to  get  a  physical  feel  for  these  four  different  solutions 
we  set  =  0  and  look  at  the  propagation  of  these  four  beams.  The  four 
incident  fields  are  listed  in  Table  3.4,  along  with  their  x  and  y 
curvatures.  From  the  table  we  see  that  these  incident  fields  all  have 
the  same  magnitudes  for  their  x  curvature  and  their  y  curvature;  the 
differences  among  the  solutions  is  in  the  signs  of  the  curvatures. 

The  beam  widths  in  the  x  direction  and  y  direction,  W  (z)  and  W  (z) , 

x  y 

are  plotted  as  a  function  of  z  in  Figures  3.2  -  3.5  for  z  values  near 
the  focal  plane.  The  following  values  were  used:  =  JI  cm., 

A  =  6328X,  f  =  10  cm.,  =  .375A,  =  .1875A  and  =  0.  Note: 

since  =  0,  8=0  and  therefore  there  is  no  need  to  use  the  rotated 
coordinate  frame.  Each  of  the  four  beams  is  quite  distinct,  except  at 
the  point  z  =  f  =  10  cm.  Here  the  four  beams  have  exactly  the  same 
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and  Wy  and  hence,  the  same  intensity  (0=0  for  all  four  beams) .  Thus 
in  practice,  if  you  have  an  astigmatic  beam  incident  upon  a  thin  lens, 
where  L>>Wq,  and  you  see  circular  isophotes  in  the  aperture  plane  and 
elliptical  isophotes  in  the  focal  plane,  you  cannot  specify  which  one 
of  the  four  beams  is  the  one  occuring  physically. 

Indeed,  since  the  second  solution  is  the  twin  of  the  first  solution, 
we  could  not  discern  between  them  even  if  the  aperture  was  smaller  (as 
long  as  the  aperture  is  circular) .  Likewise  with  the  third  and  fourth 
solutions.  If  the  aperture  is  small,  we  can  discern  between  these  two 
groups  of  solutions,  however.  This  will  be  shown  in  the  next  section. 
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a3 

B9 

a5 

B3 

B4 

B5 

a4 

°5 

a3 

°4 

“5 

Kg 

wm 

°5 

'°3 

a4 

-°5 

a3 

a4 

*5. 

K  +  “5 

a3a4/Jai  +a  5 

/ — 2 - 5 

a3a5W  °4  +a  5 

a3 

a4 

a5 

-K  *  4 

■a3a4/s /a4  +a5 

<13a5/Ja\  +a  5 

Table  3.3  All  solutions  for  ^0  and 

“4  f  0. 


Field  in  Aperture  Plane 

HHSH 

Y 

Curvature 

x  ik(a3+a4)X2  ik(a3-a4)y2 
A(p0)e  e 

1 

1 

2(ct3+a4) 

2(a3‘a4) 

...  .  -lk(“3+,14)x0  -ik(a3-!,4)5’0 
A(pQ)e  e 

-1 

-1 

2(a3+a4) 

2(a3-a4) 

,  ik(«3+a4)x0  -ik(a3-a4)y£ 
A(P0)e  e 

1 

-1 

2( “3+a4) 

_ 

2(«3-°i4) 

A,-  .  -ik(a3+a4)xj  ik(a3-a4)y2 
A(P0)e  e 

-1 

1 

2(a3+a4) 

2(a3-a4) 

Table  3.4.  Field  and  curvature  in  the  x  and  y  direction 
in  the  aperture  plane  for  solutions  that 
generate  elliptical  isophotes  in  the  focal 
piane.  (2  1  -po5/V 

A(po)  =J*wJ 


•J 


-WIDTH  IN  CM 


Z  IN  CM  *10  ! 


Figure  3.2.  Propagation  of  beam  width  in  the  x  and  y  direction  for 
solution  1  of  Table  3.3.  a-  -  .375A,  a,  =  .1875A,  a,. 
W.  -  42  cm,  f  -  10  cm  and  A  -  63288.  * 


*99-90  33-35  130-00  130-00  130.  i  3 

l  IN  CM  *10' 


o 


Figure  3.4.  Propagation  of  beam  width  in  the  x  and  y  direction  for 
solution  3  of  Table  3.3.  a.  =  .375A,  a,  =  .1873X,  a.  = 
W,..  =  JT  cm,  f  =  10  cm  and  \  *=  6328&. 
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IV.  THE  ROLE  OF  DIFFRACTION  BY  LENS  APERTURE  IN  THE  QUESTION  OF 
UNIQUENESS.  PART  II:  RADIUS  OF  BEAM _>  RADIUS  OF  LENS 

In  the  previous  section  we  have  seen  the  nonuniqueness  which 
arises  when  diffraction  from  the  lens  aperture  plays  no  role  in  the 
propagation  of  the  astigmatic  Gaussian  beam  through  a  thin  lens.  In 
the  present  section  we  are  going  to  look  at  what  happens  to  that  non¬ 
uniqueness  as  we  decrease  the  size  of  the  aperture  such  that  diffraction 
from  the  aperture  now  plays  an  important  role. 

Figure  4.1  depicts  the  physical  situation  we  are  considering  in 
this  section.  An  astigmatic  Gaussian  beam  is  incident  upon  a  thin  lens 
of  radius  L  surrounded  by  a  circular  aperture  of  variable  radius,  sL, 
where  0  <  s  £  1 . 

The  effective  radius  of  the  lens  is  then  s*L.  The  resultant 
focal  plane  intensities  generated  by  incident  fields  which  generated 
identical  focal  plane  intensities  for  the  case  Wq<<L  will  be  studied 
for  various  values  of  s. 


4.1  Calculation  of  I^(p) 

The  incident  field  we  will  be  using  in  this  section  is  the  same 
as  that  of  Section  III, 


, I Ejl 


2  2 
x  y 

~p0  /w0  ika3  2(p0//L^  _1  lka4(T2  7?)  ika52xOyO/L 

e  e  e  ^  e 


-ikculT/2  2 .  2  2  ,,  ,  2  2,  0 

3  -p0  /w0  ika3P0  ika4(x0  '  y0  }  lka52x0y0 

e  e  e  e 


(4.1) 


Figure  4.1.  Thin  Lens  Ceometry. 
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where 


u3  —  2a3/L  , 


(4.2) 


a4  =  VL  • 


(4.3) 


=  a^/L 


(4.4) 


The  aperture  plane  field  is  therefore, 


UA>  -J? t0 


-ika3L2/2  -V*  N*  ika^2  ika^x^  -  yQ2)  ik^x^ 
e  e  e  e  e 


x circ  (pQ/sL) 


(4.5) 


Looking  at  Eq.  (2.5)  and  (2.6)  we  see  that  the  focal  plane  intensity 


is  given  by 


If(p)  =  |uf(p)p  , 


Tz  | (p/xf) ! 2  , 

JTf 


(4.6) 


where 


*  co  -2iri(p/Xf)  *p0  2 

Ua(p/Xf)  -JJ  Ua(p0)e  d\ 


(4.7) 


Defining  the  spatial  frequencies  according  to 


£  -  p/Af  , 


we  see  that  in  order  to  calculate  I^(p) , 


(4.8) 
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O.  ->  CO 

ua(o  -//yp0> 


-2tt1C*P 


0  a1 

d  p. 


(A. 9) 


the  two-dimensional  spatial  Fourier  transform  of  ^(p^),  must  be  evaluated. 
Since  Eq.  (4.9),  with  Eq.  (4.5)  substituted  in,  is  not  amenable  to  ana¬ 
lytical  evaluation  we  must  evaluate  it  numerically. 

The  evaluation  of  U  (£)»  as  given  in  Eq.  (4.9),  was  performed  by 
using  a  fast  Fourier  transform  (FFT)  algorithm  to  do  the  discrete  Fourier 
transform  (DFT)  of  U  (p_) .  The  particular  subroutine  used  was  subroutine 

3  U 
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FFT  developed  by  Posey  in  1969  at  Mississippi  State  University.  This 
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subroutine  uses  the  FFT  algorithm  of  Cooley  and  Tukey  .  For  further 
details  on  the  application  of  this  subroutine  to  the  problem  we  are  con¬ 
sidering,  see  reference  21. 


4.2  Uniqueness  in  the  Phase  Retrieval  Problem 

We  are  now  ready  to  look  at  the  effect  of  the  size  of  the  aperture 
on  the  solutions  of  the  previous  section.  Throughout  this  section  we 
are  going  to  look  at  astigmatic  Gaussian  beams  of  wavelength  X  =  6328$, 
passing  through  a  thin  lens  of  radius  4  cm.  and  focal  length  10  cm.  The 
beamwidth  in  the  aperture  plane  is  taken  to  be  Wg  =  JJ  cm. 

Except  for  Figure  4.2  all  the  intensities  computed  in  this 
section  have  as  their  input  the  function  U  (p  )  sampled  at  the  points 

pg  =  (-8  cm.  +  tnA,  -  8  cm.  +  nA)  » 


m  =  0,1,2, ...,M-1 
n  =  0,1,2,. ..,N-1 


where 


and  A  =  .125  cm.,  M  =  N  =  128.  Also,  except  for  Figure  4.2,  all  the 
focal  plane  intensities  plotted  out  in  this  section  correspond  to  the 
following  focal  plane  points: 

p  =  (-2  cm  1  +  jjj,  -2  cm-1  +  jj)  \f  , 
where  M,  N  and  A  are  as  specified  above  and 

p  =  0,1,2, — ,64 
q  =  0,1,2,.. .,64 

4.2,1  Effect  of  Diffraction  on  a  Solution  with  Circular  Isophotes 
in  the  Focal  Plane 

In  this  subsection  we  compare  the  focal  plane  intensities  generated 
by  two  fields  which,  according  to  Table  3.2,  have  different  phases  but 
produce  the  same  focal  plane  intensities  if  W^<<L.  In  particular  we 
compare  the  focal  plane  intensities  generated  by 

Field  A:  =*  .75A,  =  =  0 

Field  B:  8^  =  .75A,  (5^  =  65  =  0 

for  two  different  stop  values. 

Figures  4.2a  and  4.2b  show  the  focal  plane  intensities  generated 
by  Field  A  and  Field  B,  respectively,  for  s  =  1.  In  this  case  =  J 2 
cm  and  sL  =  4  cm,  so  is  a  good  deal  less  than  sL  and  the  two  fields 
create  the  same  focal  plane  intensity  patterns.  (Note:  In  these  two 
figures  the  input  function  was  sampled  at  M  *  N  =  128  and  A  =  .0625. 
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The  focal  plane  intensities  are  plotted  over  the  region  -4Af  £  x  _<  4Af, 
-4Af  <_  y  _<  4Af . ) 

Figures  4.3a  and  4.3b  show  the  focal  plane  intensities  generated 
by  Field  A  and  Field  B,  respectively,  when  the  stop  is  reduced  to  h,- 
The  two  focal  plane  intensities  are  now  quite  different  from  each  other. 
Indeed,  Field  A  no  longer  creates  a  circularly  symmetric  intensity 
pattern. 

Clearly,  the  nonuniqueness  which  was  present  when  the  beam  was 
much  smaller  than  the  lens  aperture  is  eliminated  when  the  aperture 
plays  a  significant  role.  We  have  found  this  to  be  the  case  for  each 
of  several  examples  taken  from  Table  3.1  and  3.2  and  expect  it  to  be 
true  in  general  except,  of  course,  for  the  case  where  Field  B  is  the 
twin  of  Field  A.  It  will  be  shown  in  the  next  subsection,  however, 
that  in  practice  there  exists  a  limit  on  how  small  the  stop  can  be 
made  in  order  to  demonstrate  this  uniqueness. 

4.4.2  Effect  of  Diffraction  on  a  Solution  with  Elliptical  Isophotes 

In  this  subsection  we  compare  the  focal  plane  intensities  generated 
by  two  fields,  which,  according  to  Table  3.3,  have  different  phases  but 
produce  the  same  focal  plane  intensities  if  Wq<<L.  In  particular  we 
compare  the  focal  plane  intensities  generated  by 


Field  A: 

03  =  .375a,  =  .1875A, 

=  0 

Field  B: 

33  =  . 1875A ,  B4  =  .375A, 

s5  =  0 

for  several  different  stop  values.  Similar  results  hold  true  for  cases 
where  l4  0  and  £5  j4  0,  except  in  a  rotated  frame. 


Figure  4.3.  Focal  plane  intensity:  stop  =  \ 
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Figures  4.4a  and  4.4b  show  the  focal  plane  intensities  generated 
by  Field  A  and  Field  B,  respectively,  for  s  =  1.  In  this  case  =  JT  cm 
and  sL  =  4  cm,  so  is  a  good  deal  less  than  sL  and  the  two  fields  pro¬ 

duce  the  same  focal  plane  intensities. 

Figures  4.5,  4.6  and  4.7  each  compare  the  focal  plane  intensities 
generated  by  the  two  fields.  The  stop  values  for  the  three  sets  of 
figures  are  s  =  3/8,  %,  \  respectively.  Clearly,  due  to  the  diffraction 
from  the  aperture,  Field  A  and  Field  B  no  longer  produce  the  same  intensity 
pattern  and  the  nonuniqueness  which  was  evident  in  Figure  4.4  has  been 
eliminated. 

In  Figure  4.8  the  stop  has  been  decreased  even  further,  to  s  =  1/8. 
Here  the  two  fields  again  generate  the  same  focal  plane  intensity 
patterns,  so  there  is  some  limit  upon  how  small  we  can  make  the  stop 
and  still  demonstrate  the  uniqueness.  Physically,  what  is  happening 
in  this  case  is  that  the  aberration  and  the  stop  are  so  small  that  the 
aperture  plane  fields  for  Field  A  and  Field  B  are  almost  identical. 

It  should  be  noted  at  this  point  that  decreasing  the  sign  of 
the  stop  will  never  allow  us  to  distinguish  between  two  fields  which  are 
twins  of  each  other,  e.g.,  the  first  and  second  solution  in  Table  3  or 
the  third  and  fourth  solution  of  that  table.  More  will  be  said  in  the 
next  section  about  the  twin  solution  and  a  method  for  eliminating  it 
will  be  suggested. 


a)  a3  ■  .375A,  a,  =  .1875A,  a,  » 

b)  8,  =  . 1875A,  0.  -  . 375A ,  0^  ■ 


Figure  4.6.  Focal  plane  intensity:  stop  ■ 

a)  Ot  ■  • 375A,  a,  “  .  1 87 5 X »  a 

b)  e.  -  .  1 875 A ,  fj  -  .375A,  0 


V.  SUMMARY  AND  TWIN  SOLUTION 


In  this  section  the  "twin  solution"  ambiguity  will  be  discussed  and 
a  method  to  eliminate  it  will  be  proposed.  As  we  saw  in  Section  III,  the 
twin  solution  is  always  possible  for  the  types  of  aperture  plane  fields 
we  have  used  so  far.  In  fact,  the  twin  solution  is  always  possible  for 
any  aperture  plane  field  that  has  inversion  symmetry. 

5 . 1  Summary  and  Twin  Solution 

Phase  retrieval  problems  for  cases  where  the  incident  field  is  in 
the  form  of  an  astigmatic  Gaussian  beam  have  been  studied.  It  was  found 
that  for  cases  where  the  radius  of  the  beam  is  much  smaller  than  the 
radius  of  the  lens,  phase  retrieval  is  nonunique. 

More  precisely,  for  cases  where  the  aperture  plane  and  focal 
plane  intensities  are  circularly  symmetric,  there  exists  an  infinite 
number  of  solutions  to  the  phase  retrieval  problem  if  the  focal  plane 
intensities  have  circular  isophotes.  If  they  have  elliptical  isophotes, 
then  there  are  only  four  possible  solutions  to  the  problem,  two 
distinctly  different  solutions,  each  with  its  own  "twin". 

For  the  numerical  examples  we  did,  letting  the  lens  be  surrounded 
by  a  circular  aperture  of  significant  radius  (in  comparison  with  the 
beam  radius)  will  then  reduce  the  nonuniqueness  to  only  the  twin  solution. 

We  shall  now  discuss  the  twin  solution  further.  Let  us  now 
consider  the  incident  field 

+  +  ikW(p  ) 

UiA(p0’0)  “  lUi<P0’0)1  S  {bA) 


and  its  twin 
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-ikW(-p  ) 

UiB(P0.0)  =  IVV0)l  6  °  *  (5'2) 

The  corresponding  aperture  plane  fields  are 

UaA^0)  =  UiA(P0*°)P(V  ’  (5‘3) 

UaB^0)  =  UiB(P0’°)P(P0)  ’  (5>4) 

where  P(p^)  represents  the  pupil  function. 

It  follows  from  Eqs.  (5.1)  -  (5.4)  that 


UaA<?0>  - 

lUa<»0>l 

ikW(pQ) 

e 

(5.5) 

UaB^pO^ 

lUa<?0>l 

-ikW(-p-) 

e 

(5.6) 

where 

Ua(PQ)  =  lUi<P0*0)i  P(P0)  *  (5'7) 

Theorem 

If 

iV-^  =  ’  f°r  311  PQ  *  (5,8) 


then 


^aB(^  '[ftaA(^]*  *  for  a11 


(5.9) 


Proof: 

It  follows  from  Eqs.  (5.5)  and  (5.6)  and  the  definition  of  the 
Fourier  transform  that 


wnalfHtfiiji Hi 
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ikW(p-)  -27ri£»p-  „ 

-I'IV?o>I  *  °*  '‘V 


a„tf>  -  I * 


-ikW(-pQ)  -2niC*P0  2 


d  p, 


If  we  define  n  =  -pQ,  Eq.  (5.11)  can  be  rewritten  as 


*aB<t>  "5 |UQ(-n)|e 


-ikW(n)  2iri?-n  2 


»  '  a 


d  n  • 


Using  Eq.  (5.8)  in  Eq.  (5.12)  yields 

-ikW(n)  2Tri£*n  ,2 


«„<*>  ■  £ !».<"> I 

-  [SaA<?>  ]  *  ■ 

It  follows  from  Eq.  (2.7)  that 


d  U  , 


!fB(p)  *  IfA(p)  for  a11  p 


From  Eqs.  (5.7)  and  (5.8)  we  see  that  if 


Vv0)l  =  lui(po’0)l  * 


and 


(5.10) 

(5.11) 


(5.12) 


(5.13) 


(5.14) 


(5.15) 


P(-P0)  =  P(pQ)  ,  (5.16) 

then  | Ua (— p^q)  J  =  |ua(Pg)|  a°d  the  twin  solution  will  occur. 

Note  that  all  the  incident  fields  discussed  so  far  have  circularly 
symmetric  intensities  and  thus  satisfy  Eq.  (5.15).  Also,  the  circular 
aperture  used  satisfies  Eq.  (5.16).  As  a  result,  twin  solution  is 
always  a  possibility  in  Section  III  and  IV. 

If  the  pupil  function  is  not  of  the  type  in  Eq.  (5.16),  then  the 
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theorem  will  not  be  applicable.  Therefore,  if  we  have  a  nonsymmetric 
pupil,  the  twin  solution  might  not  occur.  In  the  next  section  we  will 
look  at  one  particular  example  where  this  is  true. 

5.2  An  Example  with  Nonsymmetric  Pupil 

A  numerical  example  of  the  same  nature  as  in  Section  IV  was 
performed  to  determine  the  effect  of  a  triangular  aperture  on  the  twin 
solution  of  the  incident  astigmatic  Gaussian  beam  with  =  0,  =  .  75A 

and  =  0,  A  =  6328$  and  =  JJ  cm.  The  aperture  used  was  an  equi¬ 
lateral  triangle  inscribed  into  the  lens  aperture  for  a  lens  of  radius 
1  cm  and  focal  length  10  cm. 

The  focal  plane  intensity  for  the  corresponding  aperture  plane 
field  is  computed  in  the  same  manner  as  in  Subsection  4.4.2  and  plotted 
out  in  Figure  5.1a  together  with  the  plot  for  its  "twin"  in  Figure  5.1b. 
It  is  evident  from  the  plot  that  for  this  particular  case,  the  use  of 
a  nonsymmetric  aperture  eliminates  the  "twin  solution"  ambiguity.  More 
studies  must  be  made  before  any  general  statement  can  be  made  about  the 
use  of  nonsymmetric  aperture  to  eliminate  the  "twin  solution". 


APPENDIX  A 


CIRCULAR  ISOPHOTES  IN  THE  FOCAL  PLANE  FOR  SOLUTIONS  IN  TABLE  3.1  AND 

TABLE  3.2 


The  intensities  in  the  focal  plane  for  an  astigmatic  Gaussian 
beam  incident  on  a  thin  lens  is  given  by  Eq.  (3.42)  with  z  =  f. 


-Zx2/W_2(f)  -Zy2/W_2(f) 

2  ^  y 

T(x,y,z=f)  =  -  «_(£)  w_(f)  '  (A,1) 

where 

W-(f)  =  WQ  1(1  +  -L-  )2T  (  )2  ,  (A. 2a) 

\l  R —  R 

xo 

w_(f)  =  WQ  1(1  +  )2  +  (  j-  )2  ,  (A. 2b) 

J  % 

Looking  at  Eq.  (A.l)  we  see  that  I^(p)  is  circularly  symmetric  if 
and  only  if 


Vf)  =  Vf) 


(A. 3) 


Substituting  Eq.  (3.36)  into  Eq.  (A. 2)  we  then  have 


Thus,  for  cases  with  =  0  (Table  3.1),  or  ■  0  (Table  3.2)  it  follows 
that 

W^f)  =  W_(f)  . 

Hence,  for  these  cases,  the  focal  plane  intensities  are  circularly  symmetric. 


S( 
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